Baryon formation and dissociation in dense hadronic and quark matter 
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We study the formation of baryons as composed of quarks and diquarks in hot and dense hadronic 
matter in a Nambu-Jona-Lasinio (NJL)-type model. We first solve the Dyson-Schwinger equation 
for the diquark propagator and then use this to solve the Dyson-Schwinger equation for the baryon 
propagator. We find that stable baryon resonances exist only in the phase of broken chiral symmetry. 
In the chirally symmetric phase, we do not find a pole in the baryon propagator. In the color- 
superconducting phase, there is a pole, but is has a large decay width. The diquark does not need 
to be stable in order to form a stable baryon, a feature typical for so-called Borromean states. 
Varying the strength of the diquark coupling constant, we also find similarities to the properties of 
an Efimov states. 



A baryon is a color-singlet bound state of three con- 
stituent quarks. Since the interaction between two 
quarks is attractive in the color-antitriplet channel, 
baryon formation can be regarded as a two-step process: 
first, two quarks combine to form a diquark with color- 
antitriplet quantum numbers, and then this diquark com- 
bines with another color-triplet quark to form a color- 
singlet bound state 

At extremely high baryonic densities and low temper- 
atures quarks form Cooper pairs in the attractive color- 
antitriplet channel, leading to the phenomenon of color 
superconductivity (loMT3j [for recent reviews, see e.g. 
Refs. 0, EB|]- Because of asymptotic freedom, the inter- 
action is weak and, just like in BCS theory, the Cooper 
pair wave function has a correlation length that exceeds 
the interparticle distance. However, as the density is low- 
ered, the interaction strength increases and the Cooper 
pair becomes more and more localized flol [l7| . Even- 
tually, Cooper pairs will form tightly bound molecular 
diquark states [18| . These may pick up another quark 
with the right color to form a color-singlet baryon. This 
is what must happen across the deconfinement transition 
into the hadronic phase. Understanding the nature of 
the transition between dense hadronic and quark matter 
is one of the scientific goals of the Compressed Baryonic 
Matter (CBM) experiment planned at the Facility for 
Antiproton and Ion Research (FAIR) (l9l |. 

In this paper we investigate the formation and dissoci- 
ation of baryons in different regions of the phase diagram 
of strongly interacting matter: the phase of broken chiral 
symmetry (hadronic phase), the phase of restored chi- 
ral symmetry (the quark-gluon plasma) above and below 
the dissociation boundary for diquarks, and the phase 
where quark matter is a color superconductor. We use 
an NJL-type model [2(1 Hl[ for two quark flavors and 
employ the following strategy. First, we compute the 
full propagator for the scalar diquark state via solving a 
Dyson-Schwinger equation. With the diquark propaga- 



tor and an additional quark propagator, we then solve a 
Dyson-Schwinger equation for the baryon propagator. 

Our approach bears some similarities to previous stud- 
ies of diquark and baryon formation (22Tl26j . These works 
also considered an NJL-type model, but they solved 
the full Faddeev equation instead of a (simpler) Dyson- 
Schwinger equation to obtain baryon states. The differ- 
ence is that in the Faddeev equation the coupling be- 
tween quark and diquark is not assumed to be local: a 
non-static quark can be exchanged between them. Our 
work is based on the cruder approximation of a local 
quark-diquark coupling. These works also considered 
the axial-vector diquark state, not only the scalar one, 
and thus were able to investigate also excited baryon 
states. On the other hand, in those works only the 
zero-temperature case was studied, while we also consider 
non-zero temperature. Moreover, we do not assume the 
diquark to be a well-defined quasi-particle in order to 
solve the Dyson-Schwinger equation (an approximation 
employed in the aforementioned works in order to solve 
the Fadeev equation). We shall see that diquarks can 
also be unstable, but still give rise to stable baryons, a 
typical feature of a Borromean state also encountered in 
atomic and nuclear physics. Varying the diquark cou- 
pling strength, we also find that our baryon has proper- 
ties which bear similarities to those of an Efimov state. 
We use natural units % = c = ks = 1; the metric tensor 
is = diag(+, -,-,-). 

The Lagrangian of the two-flavor NJL model with 
diquark-diquark interactions reads 



Cnjl = i>(ij^ 1 - ra + A7o)"0 

+G s [(4>i>) 2 + (^7 5 T^) 2 ] 

+Gd [lpij5T2 Jalpc] [V^TSTa Jai>] 



(1) 



Here, we have suppressed the color indices in the funda- 
mental representation, a = 1,2,3, and the flavor indices, 
a = u,d, in the quark spinors ip = i\) aa . The bare mass 
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matrix is mo = diag(mi°\ and the chemical poten- 

tial matrix is ft = diag(/i u , fj,^), t s (s = 1,2,3) are the 
Pauli matrices in flavor space, (J a )bc = —i e abc are the 
antisymmetric color matrices, Gs and Go are coupling 
constants for quark-antiquark and quark-quark interac- 
tions, respectively. In principle, Go can be related to 
Gs via a Fierz transformation, but we choose to keep it 
as a free parameter, allowing to explore a wider range 
of potentially interesting phenomena within our effective 
model for the strong interaction. 

In the following, we neglect the contribution from the 
isovector quark-antiquark channel, ipi^Tip = 0. We 
also decompose the scalar quark current in terms of a 
condensate part and a fluctuation, tp a ip a = cr a + S a , 
where a a = ( - Q V'c<) is the chiral condensate, and we 
work in the mean-field approximation, i.e., we neglect 
terms of order 0(S%). Similarly, we decompose the di- 
quark current as '0*75 T 2<^aV'c = (A a + S a )/(2Go) and 
drop the quadratic term in 8 a , where the diquark con- 
densate is A Q = 2Go (ipij5i~2Jaipc)- The diquark con- 
densate fluctuation can be introduced by the replacement 
A a — > A a + tp a and keeping quadratic terms in the fluc- 
tuation (p a . The above operation is equivalent to per- 
forming the Hubbard-Stratonovich transformation in the 
diquark sector. The Lagrangian JTJ) now becomes 
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Here ^ = (ip,ipc) T and ^ = (ip,tp c ) are quark spinors 
in the Nambu-Gorkov (NG) basis. The charge-conjugate 

spinors are defined by ipc — Gtp and ip c = i\> T G with 
C = ij 2 ~f°. The complex diquark fluctuation ip a has 
been decomposed in terms of its real and imaginary parts, 
<fia = {tpaR + ifai)/ V%) with color indices a = 1, 2, 3. The 
inverse fermion propagator S* -1 in the NG basis is given 
by 



PuTf + 



1I5T2 



s- 1 (p) = - 

» !b'2'-'a'-*a * fi V — p J — uv j 

(3) 

where rh = diag(m u , m^) is the quark mass matrix 
with corrections from chiral condensates, m, 



2Gs(<J u + &d) with i = u, d. The quark-quark-diquark 
vertices T ai are given by T aR = -^1^2 J a T^ G , T aI = 
-j^lbT2 J a T2 G , where (s = 1, 2, 3) are Pauli matrices 
in NG space. In the following, without loss of general- 
ity we choose the diquark condensate to be A Q = 5 a 3A3. 
Note that we only consider the scalar channel for the di- 
quark condensate, as we are only interested in the lowest 
baryon state, not the higher-lying excited ones. Includ- 
ing the axial-vector channel is straightforward, but will 
not modify our results qualitatively. Finally, we remark 
that the tadpole term f a ^ a + y*A a , which in princi- 
ple also appears in Eq. is cancelled by the term 



l Paipc'i'!§ T 2Ja'ip_+ ( PltpH5T2Ja^c at the one-loop level, 
where ipc' l J J + V-'V'c contracts and forms a quark loop in 
the NG basis. The cancellation condition is just the gap 
equation for A. 

We now add the baryon field to our Lagrangian. We 
assume the baryon to be generated by an interaction term 
between two quark and two diquark fields, 

C B = G B (pl'>pai>b'Pb 

1 -BB + ^Bf Bi tf o¥ > 



\<Pai*aT' Bi B (4) 



2G B ' 2 

Here, we decomposed ipa^a = (V'aVa) + Pen defined the 
baryonic field as B = Gb (ipafa), and neglected terms of 
order 0(/? 2 ). The baryonic fields in the NG basis are then 
denoted by B = (B, B C ) T and B = (B,B C ). The baryon- 
quark-diquark vertices are Tbr. = -^Ijvg an d ^bi = 
i^ T z G i respectively. The sum of the Lagrangians 
and Q is the starting point for our further treatment. In 
the following, for the sake of simplicity we assume exact 
isospin symmetry and we work in the chiral limit, i.e., 

a u = ad = a, thus m u = nid = m q , fi u = = fiq, and 

(0) (o) _ n 
m u = m d = (J. 

We now derive the full diquark propagator via the 

Dyson-Schwinger equation, 

A>Oo,p) = ~4£TJ -Hi.oCpo.p) , (5) 

where po = i2-KnT&xe the bosonic Matsubara frequencies 
(n = 0, ±1, ±2, . . .), i,j = R,I, and a, b = 1,2,3 are 
fundamental colors. The full propagator _D^ Q and the 
self-energy n i a only carry one index i = R,I and one 
color index a, because they are diagonal in the space of 
R, I and in color space. The self-energy has the property 
n R //, a = 5 (Ilg ± II?), where ng and Ilf depend on the 
diagonal and the off-diagonal parts of the quark propa- 
gator, respectively, and IIJ = i^nf ■ The expressions for 
ng(p ,p) and n?(p ,p) are 
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where summations over e, e', ei,ej = ±1 are implied, 
f(x) = l/(e x ' T + 1) is the Fermi-Dirac distribution, 



E k = Jk 



Ck,p+k 



Q = eE k - n, e% = v/(^) 2 +A2, and 



,k-( P +k)- 



Some simple properties of D t * are: (1) D R 3 ^ Dj g 
when A 3 ^ 0; (2) = Dj\ for any i,j = R,I; (3) 
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R, I. We also have II" = when A Q = 0. The spectral 
density for diquarks is then given by 



Pi,a(w,p) 



1 

7T 



Im Ai( w + ij ?>p) 



[ Re A,a (w + »7, p)] + [imDj (w + n?, p)] ' 



,(7) 



where we analytically continued po — > to + i 77 with real 
w and 77 a small positive number. We have similar prop- 
erties for the spectral densities as for D~^. With the 
spectral density, we can obtain the full propagator via 
the dispersion relation 
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From the Lagrangian Q the 11-component in NG 
space of the inverse baryon propagator is Sg 1 = 
-1/(2G B ) - E, where 



4 ~„ JK 



(9) 



is the 11-component of the baryon self-energy. The quark 
propagator in NG space, 5f x , is diagonal in color space. 
In the presence of a non- vanishing diquark condensate, 
S\ x = Sf-t S^. If the diquark condensate vanishes, 
Si! = Sfi = Si! and Dft a = Dj i, for any a, b. In or- 
der to evaluate E, we insert Eq. ([5]) into Eq. ©. Since 
we are interested in baryons at rest, we shall take the 
p = limit of the positive energy component of S^ 1 , 
S B ^ + (pa,p = 0) = |Tr [Sj^A+^T ] , where A£ is the 
energy projector A* = \ [1 + s (797 • p + -f M B ) /E p ], 

with = \/p 2 + M B and s = ±1. In the homogeneous 
limit, p = 0, the energy projector assumes a simple form, 
i(l + S70), which is independent of Mb- Then, 
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we obtain the spectral density as 



Ab(w,p) = - 
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RcS B ] + {u + iri,0) 



\mS B ,{u) + ir], 
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Im,S B + (u> + ir), 0) 

where we have again analytically continued po — > u + ir/. 

In our calculations for Figs. HH61 we choose the follow- 
ing parameters: G5 = 5.1 GcV~ 2 , A = 0.65 GeV (mo- 
mentum cutoff). For Figs. [T] and [BJ we vary Gd, in order 
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Figure 1: (color online) The phase diagram obtained within 
our model. For explanations see text. 



to investigate the effect of the diquark coupling constant 
on the boundaries of the diquark dissociation and the 
color-superconducting (CSC) phase and on the baryon 
formation. For Figs. [Ml we set G D = 3.11 GeV" 2 . 
This value is in the weak-coupling region, so the diquark 
is unstable in the phase of broken chiral symmetry. Nev- 
ertheless, we shall show that a quark and an unstable 
diquark can form a stable baryon in this phase. For Figs. 
3HH w e choose Gb = 10.04 GeV -1 . The baryon coupling 
constant Gb is actually the static approximation for an 
intermediate quark propagator in the Faddeev equation. 
This approximation allows us to investigate baryon prop- 
erties also at nonzero temperature and density. We fix 
Gb to obtain a baryon mass of 940 MeV in the vacuum. 

In the phase diagram of Fig. [TJ we choose four sets 
of values for temperature and quark chemical poten- 
tial, (T,fi q ) = (0.03,0.25), (0.03,0.33), (0.03,0.36), and 
(0.15,0.36), all in GeVs. They correspond to points 
A, B, C, and D. The red solid line separates the re- 
gions (indicated by xSB/xSR) where chiral symmetry is 
broken/restored; CSC denotes the color-superconducting 
phase. The blue dashed lines show the diquark dis- 
sociation boundaries for three values of the diquark 
coupling constant, Go = 3.11,3.8,4.025 (in units of 
GeV -2 ). Below a diquark dissociation line, the equation 
ReD _1 (cj, p = 0) = has a real solution ui, the so-called 
diquark pole. The corresponding regions in Fig. [T] are 
filled with light blue, green, and magenta color, respec- 
tively. These poles also exist in the CSC phases, however, 
for the sake of clarity we choose not to color the respec- 
tive regions. The CSC phases are bounded by the red 
solid line from the left and by the dash-dotted lines from 
above (from bottom to top for Gd = 3.11,3.8,4.025, re- 
spectively). Note that the diquark coupling constants 
we have chosen here are in the weak-coupling or BCS 
regime. As we increase Gd, Bose-Einstein condensation 
of diquarks could take place in the region below the dis- 
sociation lines, provided the bare quark mass is nonzero 
HM [27-32]. Note that in Ref. pjl, a vanishing decay 
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width was imposed as an additional criterion for the lo- 
cation of the dissociation boundary. 

The numerical results for the spectral densities are pre- 
sented in Figs. [2HH1 The upper panel of Fig. [2] shows the 
diquark spectral densities in the phase of broken chiral 
symmetry (point A of Fig. [1]). In the homogeneous limit 
(p = 0, red solid line), no diquark poles exist (since 
Gjj = 3.11 GcV~ 2 is too small), and the curves are 
smooth. The middle panel shows the diquark spectral 
densities in the phase of restored chiral symmetry, be- 
low the dissociation boundary, but above the CSC phase 
(point B in Fig. Q]). In the homogeneous limit, there is 
one sharp peak at uj = 0. The non-zero width of this peak 
implies that the diquark is unstable. When temperature 
grows, the diquarks dissociate, so the peak is replaced by 
a broad bump shown in the lower panel (corresponding 
to point D in Fig. [T|). In the three panels (from top to 
bottom) of Fig. Owe show pR t s, pi^ and Pi.1/2 m the CSC 
phase (point C in Fig.[T|), respectively. For p/ 3 there are 
6— function-like peaks in the range \u\ < 2A, indicating 
stable diquarks. For p^ and these peaks attain a 

small width. Also, as p increases, all peaks become wider. 
Note that the spectral densities are not odd functions of 
ui, because p q is non-zero. We see that stable diquarks 
only exist in the CSC region. Unstable diquark poles out- 
side the CSC region are actually the diquark fluctuations 
discussed in Ref. [33[- One can also see from the lower 
two panels that there are five Nambu-Goldstone (NG) 
modes which have poles at u = for zero momenta. In 
the lowest panel, there are four NG modes, i.e., the real 
and imaginary scalar fields with red and green color. In 
the middle panel, there is one NG mode for the imaginary 
scalar field with blue color. These existence of these NG 
modes is due to the validity of the following equations: 
|n J/fi>1/2 (0, 0) + ^ = and ±n /i3 (o, 0) + ^ = 0. 

In Fig. d] we show the real and imaginary parts of 
the inverse retarded Greens function for baryons (pos- 
itive energy component), again at points A,B,C, and D 
in the phase diagram of Fig. [T] In the phase of bro- 
ken chiral symmetry with m q ^ and A = (point 
A), there are no diquark condensates or resonances but 
there are stable baryon resonances: in the first panel 
(from top to bottom), we see that ReiSg, (wb,0) = 
has a solution at ljb + 3p q ps 0.94 GeV, i.e., close 
to the rest mass of the nucleoli. There is a region of 
lu b G [— 3(m, + p q ),3(m q - p q )] or M B G [—3m q , 3m q ], 
where the imaginary part ImSg^(idB,0) is very small 
(smaller than 1CV 6 GeV) in the homogeneous limit. The 
position is just inside this region, i.e., Mb < 3m q : the 
baryon weighs less than its constituents. It is therefore 
stable, although its constituents by themselves are un- 
bound, like in a Borromean state in atomic or nuclear 
physics. 

The second panel shows the case with diquark reso- 
nances but outside the CSC phase (point B). There is 
no positive energy baryon pole in this case. In the re- 
gion of higher temperatures and quark chemical poten- 
tials where chiral symmetry is restored and where there 
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Figure 2: (color online) Diquark spectral densities for dif- 
ferent values of T and p q . The upper panel corresponds to 
the point A in the phase diagram with (T, p q ) — (0.03, 0.25). 
The middle panel corresponds to the point B in the phase 
diagram with (T, p q ) — (0.03, 0.33). The lower panel is in the 
chiral symmetric phase, corresponding to the point D in the 
phase diagram with (T, p q ) = (0.15,0.36). For all panels we 
have p = pR — pi. The red solid lines are for p = 0, the 
blue dashed and brown dash-dotted lines are for p = 0.2 and 
p = 0.4, respectively. All units in GeV. 



are neither diquark condensates nor resonances (point 
D), there are also no baryon resonances and the abso- 
lute value of ImS'jj^ is very large. This case is shown 
in the third panel. In the CSC phase (point C), there 
are baryon poles but with large imaginary parts, indicat- 
ing unstable baryon resonances, as shown in the fourth 
panel. This is confirmed by a broad bump in the baryon 
spectral density in the fourth panel of Fig. 

The results for the baryon spectral density at differ- 
ent values of T and p q are presented in Fig. [S] In the 
first and second panels (from top to bottom), where 
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Figure 3: (color online) Diquark spectral densities for dif- 
ferent values of T and p q . The upper, middle and lower 
panels are pR,3, pi,z and in the CSC phase, respec- 

tively, corresponding to the point C in the phase diagram with 
(T,(iq) = (0.03,0.36). The red solid lines are for p = 0, the 
blue dashed and brown dash-dotted lines are for p — 0.2 and 
p = 0.4, respectively. All units in GeV. 



T = 0.01,0.03 GeV, we observe that the baryon spec- 
tral density hardly changes with respect to its width or 
peak position when varying the chemical potential from 
0.29 to 0.32 GeV. In the third panel with T = 0.05 
GeV the peak position shows a small increase with in- 
creasing [i q . For these larger temperatures, however, 
the width shows a dramatic increase: the curves for 
(T,fi q ) = (0.05,0.31), (0.05,0.32) GeV are not even vis- 
ible on the current scale, implying the disappearance 
of the baryon resonances. For the curves still visible 
at T = 0.05 GeV, the widths are very large indicating 
highly unstable baryon resonances. In the CSC phase 
with (T,/i 9 ) = (0.03,0.36) GeV (the fourth panel) the 
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Figure 4: (color online) The real (blue dashed) and imaginary 
(red solid) parts of the inverse propagators for baryons as 
functions of energy ui at different T and p, q . From top to 
bottom, the first panel: T — 0.03 and p q = 0.25 (point A). 
The second panel: T — 0.03 and fj, q — 0.33 (point B). The 
third panel: T = 0.15 and /j, q = 0.36 (point D). The fourth 
panel: T = 0.03 and p q = 0.36 (point C). All units in GeV. 
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Figure 5: (color online) The baryon spectral densities at dif- 
ferent values of T and \i q as functions of to + 3/j, q for p = 0. 
Four values of fi q are chosen for each panel, 0.29 GeV (red 
solid), 0.30 GeV (blue dashed), 0.31 GeV (brown dash-dotted) 
and 0.32 GeV (light blue dash-dot-dotted). In the fourth 
panel (from top to bottom) we show the result for T — 0.03 
GeV and /i, = 0.36 GeV (point C of Fig.rjJ). 
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Figure 6: (color online) The quantity lob + 3(/i 9 — m,) as a 
function of the inverse renormalized coupling 1/G r at T = 
fiq — 0. In the shaded region, i.e., below the blue curve, 
baryons are stable. The red curve inside this region is given 
by MB, P hys — 3m q , where MB,phys is the physical mass of the 
baryon. 



baryon resonance is also quite unstable, since the peak 
is very low and broad on the scale of the other panels in 
this figure. 

In Fig.[6]we vary the diquark coupling constant in order 
to investigate where the baryon is stable at T = /.i q = 0. 
We choose as x— axis the renormalized coupling G r de- 
fined in Eq. (40) of Ref. Q3. The advantage of us- 
ing G r instead oi Go is that the existence of stable 
diquark bound states is determined by the sign of G r : 
for G r > we have diquark bound states, for G r < 
they do not exist. The shaded region in Fig. [6] indi- 
cates where baryons are stable, i.e., where the imaginary 
part of the inverse baryon propagator vanishes, or where 
the spectral density may exhibit a <5-function-like peak 
(provided the real part also vanishes inside this region). 
Above the blue curve the system is in a three-quark state 
(for weak diquark coupling) or in a quark-diquark state 
(for strong diquark coupling). At moderately weak neg- 
ative G r the diquark is not stable, but, as indicated by 
the red curve, we obtain a stable baryonic bound state 
with mass Ms,phys < 3m 9 , where Ms, P hys is defined as 
the location of the peak position of the baryon spectral 
density. If we increase G r towards positive values, i.e., 
in the range where diquarks are stable, the pole energy 
of a stable baryonic bound state must lie in the range 
[— (lj d + 2fjL q ) — m q ,uj D + 2[i q + m q ] (lj d is the energy of 
the diquark at p = 0). The upper boundary of this range 
corresponds to the blue curve which is consequently given 
by ujd + 2/x g — 2m q . 

The threshold for stable baryons shown in Fig.[B]by the 
blue curve is similar to the boundary for Efimov states in 
non-relativistic cold atom physics: there, the boundary is 
proportional to — where a s is the scattering length. 
In our case, G r - a s , cf. Eq. (39) of Ref. jig. The 
curvature of the boundary in Fig. [5] indeed indicates a 
quadratic behavior as a function of G r . 



7 



100 

50 





> 

<u 

o 



-50 



-100 



T=/i f/ =0 GeV 












dj=±2m q 



-2.0-1.5-1.0-0.5 0.0 0.5 1.0 1.5 2.0 
to [GeV] 



> 

a 

o 
II 



!/5 



1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



T=fi q =0 GeV 

























4-3-2-10 1 2 3 4 
M B [GeV] 



Figure 7: The diquark spectral density (upper panel) and 
the imaginary part of the inverse baryon propagator (lower 
panel). The red solid lines are results for Gd = 5.95 GeV -2 
and the blue dashed lines are for Gd = 3.11 GeV -2 . 



The red curve for the baryon bound state was com- 
puted with a fixed coupling constant Gb- There are 
some similarities between this state and an Efimov state. 
Also there, the latter cannot form, if the two-body cou- 
pling constant is too weak, i.e., for small negative G r . 
On the other hand, for a very strong two-body coupling, 
i.e., for small positive G r , there may be a competition 
between the two-body bound state and the three-body 
bound state. There are also differences to an Efimov 
state, for instance, in Fig. [5] the baryon bound state does 
not cross the decay threshold for positive G r . We per- 
ceive this to be an artifact of a fixed quark-diquark cou- 
pling Gb in our model. In a full calculation the quark- 
diquark coupling G b should vary proportional to the in- 
verse (dressed) mass of the quark exchanged between 
quark and diquark [26[. Then, Gb will also become a 
function of the diquark coupling constant Gd- A charac- 
teristics of Efimov physics is an infinite tower of higher- 
lying excited states. In order to show that they also 
occur in our case, we would have to solve an eigenvalue 
equation for baryonic bound states. This is a subject for 
future investigations. 

In order to see the interplay between the stable di- 
quark and baryon more explicitly, we present in Fig. [7] 
the diquark spectral density and the imaginary part of 
the inverse baryon propagator for T = fj, q = GeV. For 



a strong diquark coupling (red solid line), one finds two 
components in the spectral density, a continuous compo- 
nent p c and a pole one, 

Ps(u,p) = A(p)S[u - ujpip)} - A(p)8[uj + uj p (p)}, (11) 



of 



where the amplitude is given by A(p) 
(dRen/dw)- 1 !^^), and u p (p) is the energy 
the pole with p = |p|. If the diquark coupling is 
weak (blue dashed line), only the continuous compo- 
nent remains, indicating an unstable diquark. Both 
components are taken into account in calculating 
the baryon self-energy. From the imaginary part of 
the inverse baryon propagator, one finds a region 
Mb € [— 3m q ,3m q ] where Im»S , £ 1 + = GeV in the 
weak-coupling case G D = 3.11 GcV~ 2 (blue dashed 
line), where a stable baryon can be formed. In the 
strong-coupling case Gd = 5.95 GeV~ 2 (red solid line), 
two additional bumps appear which overlap with the 
window Mb £ [— 3m q ,3m q ]. Since nonzero ImS^^ 
indicates unstable baryons, the region for stable baryons 
is reduced. This shows that the interplay between pole 
and continuum part of the diquark spectral density is 
an important ingredient in the formation of baryons. 
Neglecting the latter and taking only the pole part into 
account misses important physics (such as the formation 
of a Borromean-type stable baryon from an unstable 
diquark and a quark). 

Finally, we would like to make some comparison to 
previous works. In the Faddeev approach [22], [23|, it 
is assumed that the baryon is stable and the baryonic 
T— matrix has a separable form, which reduces the full 
Faddeev equation to the Bethe-Salpeter equation (BSE) 
for the baryonic vertex. Furthermore, for numerical sim- 
plicity it is also assumed that the diquark is stable. Thus, 
the baryon mass can be obtained via solving an eigen- 
equation (i.e., BSE) for the baryonic vertex. In this ap- 
proach, the effect of temperature was so far neglected due 
to the increase in numerical complexity. An unstable di- 
quark would also make the equation numerically hard to 
solve. Thus, so far the baryon was only treated as a stable 
bound state of a quark and a stable diquark. Since the 
baryon is stable by assumption, the properties of baryon 
resonances cannot be obtained in the Faddeev approach, 
where the baryon dissociation condition is simply real- 
ized by the condition that the baryon mass exceeds the 
sum of quark and diquark masses. However, this baryon 
dissociation condition is not correct in the CSC region 
where quarks are gapped. The correct way is to find if 
there are S— function-like peaks in the baryon spectral 
density, as done in this paper. Our static approximation 
simplifies the Faddeev equation to an RPA-type quasi- 
fermion BSE, so the baryon formation and dissociation at 
nonzero temperature and chemical potential is amenable 
to treatment. As we have shown, we have calculated the 
full baryonic spectral densities in different phases, from 
which the baryon dissociation condition is correctly ob- 
tained. 

The authors of Refs. [2J, u5fl also used the static 
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approximation in order to simplify the Faddeev equa- 
tion, but they focus on different issues. For the di- 
quark propagator, they used the proper-time regulariza- 
tion method which introduces an effective confinement, 
but the method is not applicable to nonzero temperature. 
The diquark T— matrix is approximated by a constant 
term 1/AGd plus pole terms, which is equivalent to tak- 
ing a stable diquark, while we employ the full spectral 
density of the diquark. There is some difference between 
our results and theirs. At low temperatures we also cal- 
culated the baryon mass as a function of chemical poten- 
tial: we find only a slight decrease of the baryon mass 
with chemical potential, while they obtain a significant 
decrease. The reason is that we did not include vector 
mesons and thus do not obtain large baryon number den- 
sities. Also we did not find a way of introducing confine- 
ment at nonzero temperature. We plan to look at these 
issues in a future study. In Ref. [26j, the static approxi- 
mation and a stable diquark are used. The authors con- 
sidered a three-flavor NJL model, and the baryon mass 
is found to decrease by 25% at normal nuclear matter 
density. We also plan to extend our model to the three- 
flavor case and study the properties of nuclear matter in 



the future. 

In conclusion, we used an NJL-type model to compute 
the full diquark propagator and its spectral density in dif- 
ferent regions of the phase diagram of strongly interact- 
ing matter. Baryon formation and dissociation in dense 
nuclear and quark matter is then studied via the baryon 
poles and spectral densities, incorporating the previously 
obtained diquark propagator. We find that stable baryon 
resonances with zero width are present in the phase of 
broken chiral symmetry. There are no baryon poles in 
the chirally symmetric phase. In the CSC phase, baryon 
poles exist, but they are found to be unstable due to a 
sizable width. We also pointed out that the stable baryon 
states found by us have some similarities to Borromean 
and Efimov states in atomic or nuclear physics. 
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